第 65 卷 第 3 期 天 文 学 报 Vol. 65 No.3 
2024425 H ACTA ASTRONOMICA SINICA May, 2024 


doi: 10.15940/j.cnki.0001-5245.2024.03.002 


利用 奇异 详 分 析 提 高 极 移 中 长 期 预报 精度 * 


RAT! E WI AEE 
(1 宝鸡 文理 学 院 电 子 电 气 工程 学 院 宝鸡 721016) 
(2 西安 邮电 大 学 计算 机 学 院 西安 710121 
(3 中 国 科学 院 国家 授时 中 心 西安 710600 


摘要 于 空间 大 地 观测 数据 传输 耗 时 及 处 理 过 程 复杂 , 导致 极 移 测量 值 的 获取 存在 时 延 , 无 法 满足 对 高 
精度 的 极 移 预报 值 有 重大 需求 的 应 用 领域 . 针对 极 移 复 杂 的 时 变 特性 , 提出 一 种 基于 奇异 谱 分 析 (singular 
spectrum analysis，SSA) 的 预报 方法 . 首先 用 SSA 分 离 提 取 极 移 时 序 中 的 高 频 组 分 与 低频 组 分 ; 其 次 建立 
最 小 二 乘 (least square, LS) 外 推 与 自 回 归 (AutoreGressive, AR) 模 型 对 极 移 高 频 和 低频 组 分 进行 组 合 预报 . 
结果 表明 , SSA 方法 能 够 准确 地 分 离 和 提取 极 移 低 频 和 高 频 组 分 , 对 低频 和 高 频 组 分 组 合 预报 可 以 显著 改 
善 极 移 的 中 长 期 (30-365 d) 预 报 精度 , 与 国际 地 球 自转 和 参考 系 服务 局 (International Earth Rotation and 
Reference Systems Service, IERS) 提 供 的 A 公报 中 的 极 移 预 报 值 相 比 , SSA 方 法 对 极 移 PMX 分 量 ( 本 初子 午 
线 方向 ) 和 PMY 分 量 ( 西 90° 子 午 线 方向 ) 的 中 长 期 预报 精度 改进 最 高 分 别 可 达 45.97% 和 62.44%. 研究 结果 验证 
了 SSA 方 法 对 极 移 中 长 期 预报 改进 的 有 效 性 . 
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1 引言 接 . 然而 , 由 于 数据 传输 延迟 和 复杂 的 数据 处 理 过 

地 球 定向 参数 (Earth Orientation Parameter, 程 , BOP 的 获取 存在 1-7 d 的 滞后 性 吕 ， 而 深 空 探测 
EOP) 是 实现 国际 地 球 参考 系 和 国际 天 球 参考 系 相 器 追踪 和 卫星 精密 定 轨 等 领域 对 EOP 短 期 预报 值 
互 转换 的 必需 参数 , 包括 极 移 、 世 界 时 (Universal AFEA RRP, BOP 中 长 期 预报 值 则 在 卫星 自主 
time, UT1)、 岁 差 和 章 动 . 全 球 定位 系统 (Global 。“ 导航 和 天 文 地 球 动力 学 研究 中 具有 重要 应 用 问 , 因 
Positioning System，GPS)、 甚 长 干涉 基线 测量 此 , 对 EOP 进 行 短期 和 中 长 期 预报 至 关 重 要 . 


(Very Long Baseline Interferometry, VLBI) Æ EOP 中 的 极 移 参数 定义 为 地 球 瞬 时 自转 轴 在 
载 多 普 勒 无 线 电 定位 系统 (Doppler Orbitography 地 球 本 体内 的 运动 , 极 移 在 地 极 坐 标 中 包括 PMX 
and Radio-positioning Integrated by Satellite, DO- (本 初子 午 线 方向 ) 和 PMY ( 西 90° 子 午 线 方向 ) 两 个 
RIS) 等 空间 大 地 测量 技术 可 测量 整套 的 BOP 参数 ， 分 量 , 现代 空 $ 测 地 技术 对 极 移 的 测定 精度 可 达 
建立 国际 地 球 参 考 架 与 国际 天 球 参考 架 的 直接 连 0.1 mas 量 级 . 受 大 气 、 海 洋 和 陆地 水 圈 等 多 种 
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激发 因素 的 影响 , 极 移 在 不 同时 间 尺 度 上 存在 不 
规则 变化 , 导致 极 移 的 预报 精度 随 着 预报 时 间 
的 增加 而 迅速 降低 . 极 移 提前 10 d 的 预报 精度 约 
为 4 mas, 提前 90 d 的 预报 精度 约 为 20 mas, [Al JL, 
极 移 的 时 变 特性 导致 其 预报 性 较 差 , 尤其 是 极 移 
中 长 期 预报 效果 不 理想 . 为 了 改进 极 移 的 预报 交 
R, 学 者 们 提出 多 种 预报 模型 , 例如 最 小 二 乘 (least 
square, LS) 外 推 模型 号 、 时 间 序 列 分 析 模 型 加 、 神 
经 网 络 趾 、 模 糊 推理 系统 间 、 自 协 方差 模型 站 和 联 
合 有 效 角 动量 函数 的 多 参数 极 移 模 型 90- 等 . LS 
外 推 、 自 回归 (AutoreGressive, AR) 模 型 和 自 协 方 
差 模型 易于 理解 、 建 模 简 单 ; 神经 网 络 和 模糊 推 
里 系统 属于 机 器 学 习 方 法 , 建 模 比较 复杂 , 并 
预报 效果 受 模 型 参数 影响 较 大 ; 联合 有 效 角 动量 
的 多 参数 模型 预报 效果 较 好 , 但 受 地 表 流 体 角 动 
量 预 报 数据 的 限制 , 通常 只 能 应 用 于 极 移 短期 预 
th. 为 比较 各 种 预报 模型 的 预报 性 能 , 第 1 届 EOP 预 
报 比较 项 目 (Earth Orientation Parameters Predic- 
tion Comparison Campaign, EOP PCC) 对 多 种 模 
型 的 EOP 预 报 值 进行 了 分 析 和 评估 , 比较 结果 显 
示 没 有 一 种 模型 可 以 对 所 有 EOP 参 数 和 任意 跨度 
的 EOP 预 报 精度 (平均 绝对 误差 ) 都 达到 最 优 ， ! 
方法 的 预报 精度 可 参考 Kalarus 等 423. 在 各 种 预报 
方法 中 LS 外 推 和 AR 模型 的 组 合 对 极 移 的 预报 精度 
较 高 (站, 是 应 用 比较 广泛 的 方法 , 国际 地 球 自转 和 
参考 系 服务 局 (International Earth Rotation and 
Reference Systems Service, IERS) 快 速 产品 服务 中 
心 的 A 公报 目前 就 是 利用 LS 外 推 和 AR 组 合 模 型 产 
生 极 移 预 报 值 . 

受 大 气 、 海 洋 和 陆地 水 圈 等 多 种 激发 因 
影响 , 极 移 不 仅 存 在 钱 德 勒 摆动 和 周年 摆动 , 还 存 
在 各 种 不 规则 的 短 周 期 变化 . 钱 德勤 摆动 和 周年 摆 
动 主要 影响 极 移 中 长 期 预报 , 而 短 周期 变化 主要 影 
响 极 移 短 期 预报 , 因此 , 对 于 极 移 中 长 期 预报 需要 
着 重 考虑 钱 德 勒 摆动 和 周年 摆动 . 为 提高 极 移 中 长 
期 预报 精度 , 有 学 者 采用 小 波 分 析 和 经 验 模 态 分 解 
等 时 频 分 析 技 术 对 极 移 时 序 先进 行 分 解 L37 均 , 提 
取 极 移 时 序 中 的 趋势 项 、 周 期 项 和 不 规则 项 , 然后 
再 利用 最 小 二 乘 外 推 和 自 回归 模型 对 趋势 项 、 周 
期 项 和 不 规则 项 分 别 外 推 预 测 , 研究 显示 时 序 分 
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解法 能 提高 极 移 的 预报 效果 . 奇异 谱 分 析 (singular 
spectrum analysis, SSA) 是 新 突 的 非 线 性 非 平 稳 信 
号 处 理 方法 , 它 能 够 从 小 样本 观测 信号 中 提取 有 用 
信息 , 并 且 利 用 提取 到 的 信息 对 信号 进行 周期 识 
别 、 趋 势 分析 、 消 噪 处 理 及 信和 号 重 构 05. SSA 是 经 
验 正 交 函数 的 一 种 展开 形式 , 其 优点 是 无 需 假 设 
任何 参数 模型 , 也 不 受信 号 平稳 性 条 件 的 限制 , 具 
有 识别 、 强 化 周期 振荡 信号 的 功能 , 适用 于 周期 
振荡 信号 县 加 的 信号 分 析 和 处 理 . 本 文 针对 极 移 
复杂 的 时 变 特性 , 充分 考虑 和 顾及 极 移 的 高 频 组 
分 和 低频 组 分 对 预报 的 影响 , 引入 SSA 方 法 分 离 和 
提取 高 频 组 分 和 低频 组 分 , 采用 LS 外 推 和 AR 模 型 
对 高 频 组 分 和 低频 组 分 进行 组 合 预 报 . 结果 表明 ， 
将 SSA 方 法 引入 到 极 移 预报 中 能 够 显著 降低 极 移 
的 中 长 期 (30-365 qd) 预报 的 平均 绝对 误差 . 


2 极 移 时 序 的 奇异 谱 分 析 


2.1 “奇异 谱 分 析 原 理 

设 z1,7X2,… ,ZN 是 一 组 长 度 为 N 的 一 维 观测 
时 间 序 列 , SSA 能 将 时 序 分 解 为 周期 性 、 趋 势 性 及 
拟 周 期 性 等 具有 不 同 变化 特性 的 数据 子 序 列 , 以 实 
现 对 原始 时 序 信号 识别 、 消 噪 处 理 及 信号 重 构 等 
目的 . SSA 的 主要 计算 步骤 如 下 所 述 15 

(1) #4) 32 ST ir VE FB BE. 知 窗 口 长 度 为 整数 M (1 
< M < N/2), 则 构建 的 时 清 性 矩阵 天 可 表示 为 


X — [Xi X2, ,XN-M+1 


Tı T2 ZN-M+1 


T3 TN-M+2 


M TM+Il TN 


En} yii PE 

X;= [zl za 
M +1. 

(2) 奇 异 值 分 解 . AEE RES = XXT, it BSH 

RE IE (EL, 且 将 特征 值 大 小 降序 排序 (A > Xa > 

- > Am > 0), 矩阵 对 应 的 特征 向 量 是 V1, U, 

++, Uy. #d = rank(X), Vi = XU/VM, 


抢 阵 环 的 每 行 、 每 列 为 原始 时 序 的 子 序列 ， 
,Zi IT HH, 1 <i< N- 
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二 1,2,.….,d, 则 时 灌 性 矩阵 入 可 表示 为 到 的 , 具有 极 高 的 测量 精度 , 本 文采 用 EOP 14 Cod 
序列 中 采样 间隔 为 1 d 的 极 移 数据 进行 分 析 . 利用 

X =X X -+X 2 
ERA E Aas (2) SSA 方 法 对 时 间 序列 数据 进行 分 析 时 , 重点 是 确定 
式 中 , EEX, = VApUkVE ,VA 为 矩阵 的 奇异 值 。 ”窗口 长 度 M, 通常 而 言 , 当 窗口 长 度 1 <M < N/3 
(3) 对 角 平 均 . 将 初等 矩阵 重新 转化 成 长 度 。 ”时 , 能 较 好 地 识别 周期 ML/5_M 的 周期 信号 47 对 


为 N 的 新 时 序 , 称 为 重 构 分 量 (reconstructed com- 
ponent, RC), 重 构 分 量 之 和 与 原始 时 序 相等 . 假 
WHERE Z = Xp, zj 是 经 对 角 平 均 得 到 的 重 构 分 
B,i=1,2,---, M, j = 1,2,…,N. 设 M* = 
min(M,N — M +1), K* = max(M,N —-M +1), 


EM > K*, WE; = zji, 否则 使 区) = zij, WA 
平均 公式 可 表示 为 
Zi = 
(1 = * 。 * 
和 l<i<M"*, 
i 1 zx 
| I eee M*<i<K*, 
| 三 
| 1 N-K*41 ; . 
Nil, eee a eG 
(3) 


极 移 序列 中 包含 长 期 趋势 变化 、 钱 德 勒 摆 动 、 
周年 摆动 信号 以 及 不 规则 的 短 周期 变化 信号 , 利用 
w-correlation 分 析 重 构 分 量 之 间 的 相关 性 , 将 信 


号 特征 相似 的 重 构 分 量 分 组 . 将 重 构 分 量 用 YY 个、 
YY 中 表示 , 则 重 构 分 量 之 间 的 w-correlation 可 表示 
为 19]. 
(Y®, YO) P 
a 1 < ) < M, 4 
TO yop S mS 
式 中 ， Wy | = (Y ® Y®), (YO, YO) = 


Te, wy y, uA HAM, 其 定义 为 wj 
min(k, M, N — k), HRN BY OW FAR. 
0 的 绝对 值 越 接近 于 1 BO, YO 
之 间 的 相关 性 越 强 , 应 该 把 它们 视 为 同一 信号 分 


E. 


2.2 ”分 析 结 果 
IERS 提 供 的 EOP 14 C04 序 列 是 通过 DORIS、 
GPS 和 VLBI 等 多 种 空间 大 地 观测 资料 联合 解 算 得 
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于 周期 振荡 频率 已 知 的 极 移 时 间 序 列 , 钱 德 勒 摆动 
和 周年 摆动 信号 显著 , 并 且 这 两 个 信号 是 影响 极 移 
预报 精度 的 主要 因素 , 因此 , 本 文 以 分 离 和 提取 钱 
德 勒 摆 动 和 周年 摆动 信号 为 目标 , 选择 窗口 长 度 

M = 730. 

限于 篇 幅 ， 本文 仅 以 2003-01-01 一 2014-12-31 
期 间 的 极 移 PMX 分 量 和 PMY 分 量 数据 序列 为 例 ， 
利用 SSA 方 法 对 该 时 段 序 列 进行 处 理 与 分 析 , 将 
SSA 分 解 得 到 的 前 100 阶 RC 对 应 的 奇异 值 贡献 率 
( 取 以 10 为 底 的 对 数 ) 绘 于 图 1, 奇异 值 贡 献 率 计 算 
公式 如 下 式 : 


Nk 
Nk = FE 


vA 


从 图 1 可 以 看 到 , 奇异 值 贡 献 率 较 大 的 RC 集 中 
在 前 若干 阶 , 为 更 加 直观 清楚 地 展示 极 移 PMX 和 
PMY 分 量 的 RC 奇 异 值 贡献 率 , 将 SSA 前 8 阶 RC 的 
奇异 值 贡献 率 列 于 表 1， PMX 和 PMY 分 量 的 前 8 阶 
RC 序列 分 别 如 图 2 和 3 所 示 , 这 些 RC 之 间 的 w-co- 
rrelation 分 析 结 果 如 图 4 所 示 . 从 表 1 和 图 2-3 可 以 看 

到 , PMX 和 PMY 分 量 的 1 阶 RC 的 奇异 值 贡献 率 最 
K, DI 25.43% F162.28%, 对 应 极 移 趋势 信号 ; 
2 阶 和 3 阶 RC 的 奇异 值 贡献 率 仅 次 于 1 阶 RC, 且 2 阶 
和 3 阶 RC 的 奇异 值 贡献 率 比 较 接 近 , 其 中 , PMX 
分 量 的 2 阶 和 3 阶 RC 的 奇异 值 贡 献 率 分 别 为 22.79 马 
和 21.02%, PMY 分 量 的 2 阶 和 3 阶 RC 的 奇异 值 贡献 
率 分 别 为 12.09% 和 11.29%. 从 图 4 可 以 看 到 , PMX, 
PMY 分 量 的 2 阶 和 3 阶 RC 之 间 的 相关 系数 均 大 于 
0.9. 根据 SSA 原 理 , 当 序 列 中 存在 某 一 周期 信号 时 ， 
SSA 将 其 分 解 得 到 一 对 接近 相等 的 特征 值 , 即 当 
两 个 相 邻 RC 之 间 的 相关 系数 大 于 一 ae E 
常 大 于 0.8), 则 判断 属于 同一 周期 信号 na, 因此 2 阶 
和 3 阶 RC 属 于 同一 周期 信号 , 结合 图 2- 3 可 以 判断 

这 两 个 RC 对 应 钱 德 勒 摆动 项 . PMX、PMY 分 量 


(5) 


> E 
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的 4 阶 和 5 阶 RC 的 奇异 值 贡献 率 也 比较 接近 , 它们 阶 次 越 高 , 页 献 率 越 小 , 这 些 RC 表 现 为 准 周期 性 
之 间 的 相关 系数 大 于 0.8, 同样 可 以 判断 属于 同一 荡 , 变化 规律 没有 前 5 阶 RC 那 样 明显 , 代表 亚 季 


ll 


dt 


周期 信号 , 代表 周年 摆动 项 . PMX 和 PMY 分 量 6 阶 。 ”性 准 周期 信号 . 图 2-3 中 给 出 的 各 RC 的 振荡 序列 特 
以 后 的 RC 的 奇异 值 页 献 率 相对 较 小 , EAA F, 征 与 表 1 中 列 出 的 奇异 值 页 献 率 对 应 得 很 好 . 
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Fig.1 Distribution of the eigenvalue contribution ratio of the RC 
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Fig. 2 


图 2 极 移 PMX 分 量 2003-01-01 一 2014-12-31 数 据 序 列 的 前 8 阶 RC 


First eight order RC of the PMX component time-series of polar motion between Jan 1, 2003 and Dec 31, 2014 
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图 3 极 移 PMY 分 量 2003-01-01 一 2014-12-31 数 据 序列 的 前 8 阶 RC 


Fig.3 First eight order RC of the PMY component time-series of polar motion between Jan 1, 2003 and Dec 31, 2014 
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SSA 前 8 阶 RC 的 w-correlation 
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Fig.4  w-correlations between the first eight order RC with SSA 
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#1 SSA 前 8 阶 奇异 值 贡献 率 如 图 5-6 所 示 . 从 图 5-6 可 以 看 到 , OR AST RCA IM 
Table 1 Contribution ratio of the first eight 后 构成 的 极 移 低 频 组 分 与 原始 极 移 序 列 符合 得 很 
a 好 , 重 构 低频 组 分 能 够 表现 原始 极 移 序列 的 整体 变 
a a 化 趋势 ; 剩余 RC 之 和 构成 的 极 移 高 频 组 分 表现 为 
1 2 62.28 不 规则 的 非 线 性 变化 , 适合 于 利用 时 间 序 列 分 析 模 
2 22.79 12.09 型 对 其 进行 预测 . 为 检测 极 移 低 频 组 分 和 原始 极 移 
3 21.02 11.29 序列 中 的 周期 振荡 信号 , 利用 快速 傅 里 叶 变 换 (Fast 
4 7.57 4.18 Fourier Transform, FFT) 对 极 移 低频 组 分 和 原始 
5 6.83 3.65 极 移 序列 进行 频谱 分 析 , 结果 也 在 图 5-6 中 给 出 . 从 
igh Lös 图 5-6 可 以 看 到 , 原始 极 移 序 列 和 低频 组 分 中 都 包 
~~ e 含 周 期 为 438 d 的 钱 德 勒 摆动 和 周期 为 365 d 的 周年 

| | 摆动 信号 , 表明 低频 组 分 中 的 主 周期 振荡 信号 与 原 


0:91 0.38 始 序列 相同 ; 从 频谱 图 中 还 可 以 发 现 , 除了 极 移 序 

列 中 的 钱 德 勒 摆动 和 周年 摆动 信号 , 极 移 序列 中 还 

考虑 到 前 5 阶 RC 的 奇异 值 贡献 率 较 大 且 属 于 ” 含有 振幅 较 小 、 周 期 约 为 500 d 的 年 际 信号 (9-20. 

低频 信号 , 故 将 前 5 阶 RC 线 性 又 加 后 构成 极 移 低 频 ERANA REN, SSA 方 法 能 够 很 好 地 分 离 和 提 
组 分 , 而 将 剩余 RC 线 性 受 加 后 构成 极 移 高 频 组 分 ， ” 取 极 移 序列 中 的 低频 组 分 和 高 频 组 分 . 
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(b) Low- and high-frequency components 
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图 5 极 移 PMX 分 量 的 观测 时 序 、 高 频 分量 、 低 频 分 量 及 频谱 . (a) 观 测 时 序 ; (b) 高 频 和 低频 分 量 ; (c) 观 测 时 序 的 频谱 ; (qd) 低 频 分 量 的 频谱 . 


Fig.5 Observed time-series, high- and low-frequency components of the X component of polar motion and their frequency 
spectrum. (a) Observed time-series; (b) High- and low-frequency components; (c) Frequency spectrum of observed time-series; 


(d) Frequency spectrum of low-frequency components. 
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图 6 极 移 PMY 分 量 的 观测 时 序 、 高 频 分 量 、 低 频 分 量 及 频谱 . (a) 观 测 时 序 ; (b) 高 频 和 低频 分 量 ; (c) 观 测 时 序 的 频谱 ; (d) 低 频 分 量 的 频谱 


Fig.6 Observed time-series, high- and low-frequency components of the Y component of polar motion and their frequency 
spectrum. (a) Observed time-series; (b) High- and low-frequency components; (c) Frequency spectrum of observed time-series; 


(d) Frequency spectrum of low-frequency components. 


3 ”预报 模型 之 和 . 图 7 给 出 了 基于 SSA 的 极 移 预 报 流程 
与 采用 SSA 直 接 对 极 移 进行 预报 的 方法 不 同 ，。 3.1 ， 最 小 二 乘 外 推 模型 
PERE E E E 对 于 极 移 低频 组 分 , LS 外 推 模型 包含 线性 趋势 


利用 SSA 能 够 分 离 和 识别 复杂 信号 的 周期 分 量 , 特 项 、 钱 德 勒 摆动 项 和 周年 摆动 项 , 先 验 模型 的 数学 
别 是 极 大 地 改善 周期 分 量 的 分 辩 率 的 优良 特性 , SARK 

建 LS+AR 组 合 模型 对 低频 组 分 和 高 频 组 分 分 别 预 

测 , 进一步 提高 极 移 预报 精度 . 基于 SSA 的 极 移 预 f(t) = a + bt + ecos(wit) + disin(wit) + 
报 流程 如 下 : 首先 利用 SSA 分 离 和 提取 极 移 序 列 中 
的 低频 组 分 (包括 钱 德 勒 摆 动 项 、 周 年 摆动 项 和 长 
期 趋势 项 ) 和 高 频 组 分 , 然后 针对 极 移 序 列 中 的 低 AP, t=1,2, ,N, N 表 示 极 移 观 测序 列 的 长 度 ， 
频 组 分 建立 一 个 LS 外 推 模型 , 具体 方法 为 先 利 用 一 wi 和 ws 分 别 表 示 钱 德 勒 摆动 角 频 率 和 周年 摆动 角 
段 低频 组 分 序列 进行 最 小 二 乘 拟 合 求解 先 验 模型 频率 ; 和 1 表示 极 移 的 趋势 参数 ; cl 和 qi 表示 钱 德 
参数 , 得 到 低频 组 分 的 外 推 值 及 最 小 二 乘 残 差 , 然 勒 摆动 参数 ; c, 和 ds 表示 周年 摆动 参数 ; a b cn 
后 运用 AR 模型 对 高 频 组 分 及 最 小 二 乘 残 差 之 和 建 di1、c2 和 ds 这 6 个 未 知 参数 可 以 通过 最 小 二 乘法 进 
模 预 测 , 极 移 最 终 预 报 值 为 LS 外 推 值 和 AR 预测 值 行 求解 . 


C2 cos(@wat) + də sin(@wet) 3 (6) 
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Fig.7 Flowchart of polar motion prediction based on the SSA 


3.2 BEIR 


假设 平稳 随机 时 序 为 {wt = 1,2,.… , N}, 自 


回归 模型 描述 了 当前 时 刻 的 随机 变量 值 与 前 p 个 
随机 变量 值 的 线性 组 合 及 误差 项 之 间 的 函数 关系 ， 
2 阶 自 回归 模型 的 数学 表达 式 为 D 
p 
= 5 Pilti + Et, (7) 
i=l 
式 中 , 6 表示 高 斯 白 噪 声 ; y; 表 示 自 回归 模型 的 系 


数 , i=1,2,--- ,p, 本 文通 过 求解 尤 尔 -沃克 方程 确 
定 模 型 系数 . 
自 回 归 模 型 的 阶 数 pz 关乎 AR 模型 的 泛 化 能 


本 文 利 用 赤 池 信息 论 准 则 (Akaike’s Information 
Criterion, AIC) 优 化 选取 自 回归 模型 的 阶 数 P3: 
2p 

入， (8) 
式 中 , op 代表 自 回归 模型 的 拟 合 残 差 的 标准 差 , 当 p 
二 1,2,--- , VN, AIC(p) 取 得 最 小 值 时 所 对 应 的 p 即 
为 自 回归 模型 的 最 佳 阶 数 . 


A ”模型 验证 
4.1 ”精度 评估 指标 
利用 预报 绝对 误差 (absolute error, AE) iik 
下 均 绝对 误差 (mean absolute error, MAE) 对 预报 
精度 进行 分 析 评估 , 预报 AE、 预 报 MAE 的 计算 方 
法 分 别 表示 为 3 


AIC(p) = Ino? 十 


a 


AE = |P® — 0, (9) 


1 n 
MAE, = z 2, AEn c, (10) 
式 中 , 代表 预报 跨度 ; ”代表 预报 期 数 ; P, OW 


分 别 代表 c 点 的 预报 值 、 观 测 值 . 


与 LS 十 AR 模 型 的 比较 

选择 2003-01-01 一 2019-10-25 期 间 的 极 移 数据 
进行 预报 试验 , 其 中 , 2015-01-01 一 2019-10-25 为 预 
报 期 , 每 隔 7 d 预 报 一 期 ， 基础 建 模 数 据 的 长 度 为 
12 yr, 每 期 预报 后 建 模 数据 也 随 之 更 新 , 一 共 进 行 
了 200 期 预报 . 
为 验证 SSA 方 法 对 LS 外 推 和 AR 组 合 模 型 的 极 
移 预 报 改善 情况 , 分 别 利 用 LS+AR 模 型 和 SSA+LS 
十 AR 模型 对 极 移 进 行 1-365 d 的 预报 , 图 8 为 两 种 模 
型 在 2015-01-01 一 2019-10-25 期 间 的 预报 AE. 从 图 8 
可 以 看 到 , 两 种 方法 的 最 大 预报 AE 均 出 现在 2017 
年 , 且 2017 年 的 预报 AE 明 显 大 于 其 他 年 , 这 可 能 与 
2017 年 发 生 的 La Nina (拉尼娜 ) 事 件 有 关 ; LS+AR 
模型 对 极 移 PMX 和 PMY 分 量 的 最 大 预报 AE 分 别 
超过 60 mas 和 90 mas, 而 SSA 十 LS 十 AR 模 型 对 极 移 
PMX 分 量 的 最 大 预报 AE 小 于 50 mas, 对 极 移 PMY 
分 量 的 最 大 预报 AE 小 于 80 mas，SSAT+LS+AR 
模型 的 预报 AE 明 显 小 于 LS+AR 模 型 ; LS+AR 模 
a 含有 周期 性 效应 , 说 明 受 极 移 高 频 
分 量 的 影响 , 直接 对 极 移 序列 进行 最 小 二 乘 拟 合 
无 法 充分 反映 钱 德 勒 项 和 周年 项 的 变化 特征 , 从 
而 在 预报 AE 中 出 现 残余 的 周期 项 拟 合 和 外 推 误差 ， 
而 SSA+LS+TAR 方 利用 不 同 模型 对 极 移 低 频 分 量 

和 高 频 分 量 分 别 预测 能 够 降低 预报 AE. 


4.2 
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8 ”LS 二 AR 和 SSA 十 LS 十 AR 模 型 的 预报 AE 


Fig.8 AE of the prediction of the LS+AR and SSA+LS+AR models 


为 定量 分 析 SSA 十 LS 二 AR 模型 相对 于 LS 二 AR 
模型 的 极 移 预 报 改善 效果 , 表 2 给 出 了 LS+AR 模 
型 和 SSA+LS+AR 模 型 在 不 同 预报 跨度 上 的 MAE. 
从 表 2 可 以 发 现 , SSA+TLS+AR 模 型 相对 于 LS 上 +AR 
模型 在 各 个 跨度 上 预报 精度 均 有 一 定 程度 的 提高 ， 
其 中 , 对 于 1-10 d 跨 度 的 极 移 超 短期 预报 , 精度 提 
高 在 10% 以 内 , M30 d 跨 度 开 始 , SSA+LS+AR HEAL 
的 预报 精度 提高 比较 明显 , 对 极 移 PMX 和 PMY 分 
量 的 预报 精度 最 大 分 别提 高 32.89% 和 23.14%, H. 
PMX 分 量 的 预报 精度 提高 维持 在 20% 以 上 , PMY 
分 量 的 预报 精度 提高 维持 在 15% 以 上 . SSA+LS+ 
AR 模 型 对 于 极 移 中 长 期 预报 精度 改善 更 为 显著 ， 
说 明 钱 德 勒 摆动 和 周年 摆动 信号 主要 影响 极 移 中 
长 期 预报 , 验证 了 SSA 分 离 和 提取 钱 德 勒 摆动 和 周 
年 摆动 信号 的 有 效 性 . 


ae 
ah 


i 


4.3 与 IERS A 公报 的 比较 

为 进一步 验证 SSA 模 型 对 极 移 预 报 的 改善 效 
R, 将 SSA+TLS+TAR 模 型 的 预报 结果 与 IERS AZ 
报 2019 年 第 1-52 周 提供 的 极 移 预 报 结果 进行 对 比 
分 析 , 预报 时 段 为 2019-01-01 一 2020-12-25. 以 IERS 
EOP 14 C04 序 列 作为 参考 , 分 别 计 算 本 文 预报 结 
果 和 IERS A 公报 预报 结果 的 平均 绝对 误差 , 结果 
如 图 9 所 示 , 图 9 中 蓝 色 实 线 表示 IERS A 公报 的 平 
均 绝 对 误差 , 蓝 色 虚线 表示 SSA+LS+AR 模 型 的 平 
均 绝对 误差 , 红色 实 线 表 示 SSA+LS+AR 模 型 相对 
于 IERS A 公报 的 平均 绝对 误差 改善 率 . 从 图 9 可 以 
看 到 , 对 于 1-30 d 的 短期 预报 , 本 文 方法 的 预报 精 
度 低 于 IERS A 公报 的 结果 ; 对 于 30 d 以 上 的 中 长 期 
预报 , IERS A 公报 的 极 移 预 报 误 差 随 着 预报 跨度 
的 增加 而 快速 增 大 , 误差 发 散 很 快 , 而 本 文 方法 的 
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预报 结果 虽然 在 预报 前 期 不 如 IERS A 公报 的 结果 ， H, 极 移 PMY 分 量 的 预报 效果 优 于 PMX 分 量 , 其 主 
HE30 d 以 后 的 预报 中 精度 较 高 , 预报 精度 显著 优 要 原因 是 极 移 PMY 分 量 的 数据 规律 性 更 强 , 此 外 ， 
FIERS A 公报 的 结果 , 对 极 移 PMX 和 PMY 分 量 的 极 移 预 报 结果 和 预报 时 段 也 存在 着 一 定 的 关系 . 


预报 精度 最 大 分 别提 高 45.97% 和 62.44%. 在 本 文 


Wl 


表 2 LS+AR 和 SSA+LS+AR 模 型 在 不 同 预报 跨度 上 的 MAE (mas) 
Table 2 Statistics of the MAE (in unit of mas) for the LS 十 AR and SSA 十 LS 十 AR models under 


different prediction lengths 


LS+AR SSA+LS+AR Improved percentage 
Prediction length/d 


PMX PMY PMX PMY PMX PMY 

1 0.26 0.21 0.26 0.21 0.00% 0.00% 
5 1.64 1.18 1.58 1.13 3.66% 4.24% 
10 3.23 2.15 3.03 1.99 6.19% 7.44% 
30 8.17 6.29 6.53 5.38 20.07% 14.47% 
60 15.73 14.00 11.15 10.76 29.07% 23.14% 
90 20.83 20.86 13.98 16.09 32.89% 22.87% 
120 22.77 24.06 15.74 18.90 30.87% 21.45% 
180 24.59 25.75 17.07 20.40 30.58% 20.78% 
240 27.91 27.94 20.75 22.04 25.65% 21.12% 
300 26.53 28.95 20.61 23.70 22.31% 18.13% 
365 25.83 28.75 19.03 24.17 26.34% 15.93% 
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图 9 SSA 十 LS 十 AR 模 型 和 IERS A 公报 的 预报 MAE 


Fig.9 MAE of the prediction of the SSA+LS+AR model and IERS Bulletin A 
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Improvement in the Medium- and Long-term Prediction 
Accuracy of Polar Motion Using Singular Spectrum Analysis 
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Asstrract Polar motion cannot be determined in real time owing to the delay caused by data transfer 
and heavy computation procedures. Polar motion predictions are therefore required for many real-time 
applications and geodynamics study. Polar motion with respect to the axis of the terrestrial reference 
system is not constant in time but changes due to external forces and internal processes. This paper 
proposes a hybrid method to improve polar motion prediction based on singular spectrum analysis (SSA). 
In this method, the SSA is employed to separate and extract the low- and high-frequency components 
of polar motion. Next, the least square (LS) extrapolation and autoregressive (AR) models are applied 
to model and predict the extracted low- and high-frequency components, respectively. The subsequent 
predictions of the low- and high-frequency components are summed to obtain the predicted values of 
polar motion. The singular spectrum analysis of the polar motion time-series shows that the low- and 
high-frequency components can be accurately extracted from the original time-series of polar motion by 
the SSA technology. Further, the prediction results illustrate that the proposed SSA-based combined 
model can substantially improve the medium- and long-term predictions of polar motion out 30 days 
in the future. In comparison with the polar motion predictions from the Bulletin A published by the 
International Earth Rotation and Reference Systems Service (IERS), the highest improvements found for 
PMX in the prime meridian direction and PMY in the west 90° meridian direction are 45.97% and 62.44% 
up to 365 days in the future, respectively. It is concluded that the SSA is a potential method to enhance 
medium- and long-term predictions of polar motion. 


Key words astrometry: reference systems, Earth rotation: polar motion, methods: data analysis 


24-12 


